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A review is given of the QCD corrections to charm quark production 
in deep inelastic electron-proton scattering. An outline of the computa- 
| tion of the virtual photon-parton subprocesses, from which one obtains the 

iy) . heavy quark coefficient functions, is given. The dominant production mech- 

' anisms are discussed. Further we show that the asymptotic heavy quark 

coefficient functions, computed in the limit Q 2 3> m 2 , can be derived using 
the operator product expansion technique. Further we present the various 
' schemes proposed in the literature to describe the charm component of the 

structure function and compare them with the most recent data from the 
' experiments carried out at HERA. 

a: 
a: 

3^ 1 1. Introduction 

The study of heavy quarks and their decay and production mechanisms 
provides us with important insights in the standard model of the electroweak 
$_i ' and strong interactions. In this model the heavy flavours are given by the 

charm (c), bottom (b) and the top (t) quarks. Because of the confinement 
property in the theory of the strong interactions, given by quantum chromo- 
dynamics (QCD), these quarks can only be observed via the heavy baryons 
and mesons in which they are confined. Since mesons are more copiously 
produced than baryons we concentrate on the former particles. Here one 
can distinguish between the following type of mesons. 

1 open heavy quark mesons 

In this case the heavy quark is accompanied by a lighter quark mostly 
represented by the up (u), down (d) and strange (s) quarks. Examples 
are: 

* lectures presented at the XXXVIIth Cracow School of Theoretical Physics, 30th May 
- 10th June 1997, Zakopane, Poland. 

(1) 



2 



ZAKOP PRINTED ON FEBRUARY 1, 2008 



a open charm e.g. D u = uc, = dc, D s = sc. 
b open bottom e.g. B u = ub, Bd = db, B s = sb, B c = cb 
including their anti- mesons. 
2 hidden heavy quark mesons 

Here the heavy quark is always bound to its anti-quark. Examples are: 

a hidden charm e.g. J/ip = cc 

b hidden bottom e.g. T = bb 

including the higher excitations. 

In these lectures we limit ourselves to open heavy quark production. The 
mesons are observed via their electroweak decays. An example is the decay 
of the meson D u given by D u — > K~ + /i + + which proceeds via the 
partonic reaction c — > s + fi + + where // + , v„ emerge from the virtual 
VF + -boson which is exchanged between the quarks and the lepton pair. 
There are numerous reasons why open heavy quark mesons are of interest. 
Here we will mention some of them (for a review see [jl]]). 

a The observation of rare decays in particular those of the B-mesons. 
With rare we mean decays that are suppressed in the standard model. 

b The measurement of the Cabibbo-Kobayashi-Maskawa matrix elements 
denoted by Vij = u,d - ■ ■ t) via the production or decay of the heavy 
quark mesons. A classical example is the measurement of the quantity 
V cs = cos 9 C in the process e+p — > v e +D+' X' ('X' is any inclusive final 
state). For this reaction the dominant partonic subprocess is given by 
e + s — > v e + c which proceeds via the exchange of a W-boson. 

c The study of DD and Si5-mixing in connection with CP- violation in 
the B-system. 

d The study of production mechanisms of heavy quarks provides us with 
new tests of QCD. It also enables us to measure some of the parton 
densities (see below) in kinematical regions which are not accessible in 
other types of processes. An example is the gluon density which can 
be measured via charm production at the HERA-collider. 

In these lectures we will only study those production mechanisms of heavy 
quarks where the methods of perturbative QCD can be used. Reactions 
where only nonperturbative methods can be applied like e.g., diffraction 
will not be treated here. 

In perturbative QCD physical quantities can be expanded in the strong 
coupling constant a s (/x 2 ), where the scale fj, has to be large, provided we are 
dealing with so called hard processes. The latter are described by the parton 
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model including higher order QCD corrections. The hard processes are 
characterized by the property that all kinematical invariants, on which the 
quantities depend, become asymptotic whereas their mutual ratios are fixed. 
These ratios should neither become zero nor infinite. In the case of heavy 
quark production the above condition implies that the heavy quark mass 
m should become asymptotic too. Unfortunately these types of processes 
cannot be completely described by perturbative methods. All cross sections 
also contain nonperturbative parts. The latter are represented by the parton 
densities and the fragmentation functions whose properties will be discussed 
below. 

Let us first enumerate the various reactions, including the basic partonic 
subprocesses, in which the heavy mesons are produced in a semi inclusive 
way. 

1 hadron-hadron collisions 
Example: P + P^B + B+'X' 

Lowest order partonic subprocesses: g + g — > b + b, q + q — > b + b 

2 lepton-hadron collisions 

Example: e + P^D + D+'X' 

Lowest order partonic subprocess: 7* + g — > c + c 

3 photon-hadron collisions 

Example: j + P^D + D+'X' 

Lowest order partonic subprocess: 7 + g — ► c + c 

4 photon-photon collisions 

Example: -j + j^D + D+'X' 

Lowest order partonic subprocess: 7 + 7 — ► c + c 

5 electron-positron collisions 

Example: e+ + e~ -> B + B +' X' 

Lowest order partonic subprocess: e + + e~ — > b + b 

Notice that in reaction 4 one or both photons can be virtual (indicated by 
*). In these lectures we only concentrate on reaction 2 where the charm 
(anti-) quark is produced in deep inelastic lepton-hadron scattering. We 
show this process in more detail in Fig. [I], where the momenta of the 
virtual photon and the proton are denoted by q and p respectively. Notice 
that q is spacelike (i.e. q 2 = — Q 2 < 0). Further we assume that Q 2 is small 
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Fig. 1. Production of the mesons D and D in deep inelastic electron- proton scat- 
tering via the photon-gluon fusion mechanism. 



enough so that the neutral current process with the exchange of a Z-boson 
is suppressed. The process becomes deep inelastic when Q 2 and p ■ q are 
large whereas the scaling variable x = Q 2 /2p ■ q is kept fixed (i.e. < x < 1 
but x 7^ and x ^ 1). When these kinematical conditions are satisfied then 
according to the parton model the proton can be viewed as a collection of free 
quarks and gluons (called partons). Each of them can be involved in the hard 
scattering without being influenced by the other partons which are called 
spectators. The latter produce the hadronic final state ' X' . Integration over 
the momenta and summation over the quantum numbers of the spectators 
provides us with the probability that the gluon in Fig. [l] emerges from the 
proton with the fraction z of the proton's momentum. This probability is 
given by the parton density indicated by fg{z). This function can only 
be computed using nonperturbative methods in QCD. Unfortunately these 
methods are not available yet so that one either has to resort to models for 
fg{z) or one has to obtain this density from the data. The same holds for 
the quark parton densities denoted by fq(z). The incoming gluon in Fig. |l| 
is involved in a hard scattering with the photon so that a charm anti-charm 
quark pair is produced in the final state. The charm and anti-charm quark 
pick up from the vacuum a light anti-quark and quark respectively. The 
fusion of the charm quark with the light anti-quark leads to the production of 
a D-meson. This process is described by the fragmentation function T>®(z) 
where z denotes the fraction of the momentum of the charm quark carried 
away by the D-meson. An analogous description holds for the production 
of the anti D-meson coming from the anti-charm and the light quark. The 
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Fig. 2. Kinematics of charm production in deep inelastic electron-proton scattering 



light (anti-) quarks, which are called spectators, have very low momenta so 
that they cannot change the magnitude as well as the direction of the three- 
momentum of the charm quark. Therefore the momentum of the charm 
quark can be reconstructed from the momentum of the D-meson. Like 
the parton densities the fragmentation functions can be only determined 
by using nonperturbative methods. However the latter are not at that 
stage that one can compute these functions. Therefore they can either 
be determined by models or they have to be extracted from the data. In 
spite of the poor knowledge about these functions one can assume that the 
parton densities f^(z) (k = q,g) and the fragmentation functions T)^(z) 
(a = D,B;k = c,b) are universal and process independent. The former 
only depend on the type of parton (q or g) and the type of hadron from 
which the parton emerges. The same holds for the fragmentation function 
which only depends on the heavy quark and the type of meson into which 
the former fragments. This universality remains unaltered even after QCD 
corrections to the hard processes have been included. It means that when 
these phenomenological functions are obtained from a certain process A one 
can use them as input for process B in order to make absolute predictions 
for the cross section of the latter reaction. This statement is only correct if 
the QCD corrections are carried out up to the same order in both processes. 



Charm quark production in deep inelastic electron-proton scattering pro- 



2. Electroproduction of heavy quarks 




(1) 
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Here ' X' denotes any inclusive hadronic state which means that we have 
summed over all quantum numbers and integrated over all momenta of 
the hadrons belonging to this state. Further we consider neutral current 
processes only so that the intermediate vector boson V either stands for the 
photon or for the Z-boson. The momentum transfer in the above process is 
spacelike so that we have 

q 2 = -Q 2 < with q = h-l 2 . (2) 

In the subsequent part of these lectures we will assume that Q 2 < Mf so 
that the above reaction is dominated by the one-photon exchange mech- 
anism. The computation of the cross section of reaction (1) involves a 
five-dimensional integral. The kinematical integration variables are J| 

S = [p + q) , cos $ = -=, r — — - , (3) 

\h x h \ ■ \v x vi ) 



T\ = T — m 2 = (p — pi) 2 — m 2 , U\ = U — m 2 = (q — pi) 2 — m 2 . (4) 
Further we define the scaling variables 

Q 2 p-q 

x = 7i — » y = — r • 5 ) 

2p-q p-h 

Integration over $ yields the following cross section 

d 4 a(T 1 ,U 1 ,x,y) _ 
dx dy dT\ dU\ 

a r-,.. d 2 a L (T 1 ,U 1 ,x) 2 d 2 a T (T 1 ,Ui, x" 



2nxy 



where ox and ut stand for the longitudinal and transverse photon cross 
sections respectively. In addition to the variables indicated above a and 
(k = T,L) also depend on Q 2 Eq. (2) and the heavy quark mass m. The 
factor in front of d 2 ax in Eq. (6) reflects the vector nature of the photon. 
From d 2 a k /dT 1 dUi in Eq. (6) one can compute the double differential cross 
sections [|| 

d 2 a k d 2 a k 



dY d | pi± | dXj d | pi± | 
The rapidity Y and the momentum fraction Xf are defined by 



(8) 
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In the above equations the longitudinal and the transverse momenta are 
indicated by the indices I and _L respectively. In the subsequent part of 
these lectures we are only interested in the integrated cross sections given 
by 

a k (x,Q ,m ) = J dT x J dUi — — . (9) 

The total longitudinal and transverse cross sections are related to the cor- 
responding deep inelastic structure functions Fk (k = T, L) as follows 

a L (x,Q 2 ,m 2 ) = ^^-F L (x,Q 2 ,m 2 ) 

a T (x,Q 2 ,m 2 ) = — i^Q 2 ,™ 2 ). (10) 
p-q 

Instead of Fi the experimentalists measure the structure function 

F 2 (x,Q 2 ,m 2 ) =2xF 1 {x,Q 2 ,m 2 )+F L {x,Q 2 ,m 2 ). (11) 

With the above definitions the deep inelastic scattering cross section can be 
written as 

( ' 2a 2 ™\ P \{l + (l-y) 2 }F 2 (x,Q 2 ,m 2 )-y 2 F L (x,Q 2 ,rn 2 )\ . (12) 



dx dy Q 2 



Here v S e p denotes the centre of mass energy of incoming electron-proton 
system. When the charm component of the structure function is under study 
we will add an index c to the structure functions so F k is replaced by i*fc iC - 
Before we continue it is important to emphasize that perturbative QCD 
can only predict the (Revolution of the structure functions but not their 
x-dependence. The latter is partially determined by the parton densities 
which, as we have mentioned before, are of a nonperturbative origin. 
For charm production in perturbative QCD one can distinguish between 
two different production mechanisms. Let us denote the proton state in the 
Fock space by 

| P) = a± | uud) + <i2 | uu, uud) + 03 | dd, uud) + 04 | ss, uud) 
t 

+ a H I HH,uud) . (13) 

H=c 

The two possible production mechanisms are: 
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Fig. 3. Feynman diagrams for the lowest-order photon-gluon fusion process con- 
tributing to the coefficient functions . 



A Intrinsic heavy quark production [Q] 

Here we have an ^ for at least one H with H = c,b, t. In the case 
of intrinsic charm we have a c 7^ and a& = 0, at = 0. The main 
production mechanism is given by the flavour excitation process 

7* + c^c, (14) 

which is a zeroth order process in the strong coupling constant a s . 
B Extrinsic heavy quark production^] 

Here we have an = for H = c, b, t. The dominant production mech- 
anism is now given by the photon-gluon fusion process. On the Born 
level (first order in a s ) it is given by (see Fig. ||) 

7* + 3^c + c. (15) 

Notice that in the case of extrinsic charm production only light quarks 
(u,d,s) and the gluon can appear in the initial state of the partonic processes 
since the probability that a charm quark emerges from the proton is zero. In 
these lectures we will limit ourselves to extrinsic charm production because 
the recent experiments at HERA ||,(7j indicate that the data favour this 
production mechanism so that = 0. 

3. Exact heavy quark coefficient functions up to order 

The calculation of the partonic cross sections denoted by cr^fc is straight- 
forward (see ||). After using the same relations as presented for the struc- 
ture functions in Eq. (12) we obtain the heavy quark coefficient functions 
represented by Hi & (i = 2, L; k = q, g). The lowest order contribution to the 
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Fig. 4. Virtual gluon corrections to the proces 7* + g — > c + c contributing to the 

(2") 

coefficient functions i?<„. 



structure functions corresponding to the graphs in Fig. (^) can be written 
as 

F^(x,Q 2 ,m 2 ) = e 2 c J^^f^x/z)HQ(z,Q 2 ,m 2 ,a s ) 

= (16) 

with z max = Q 2 /(Q 2 + 4m 2 ) and e c denotes the charge of the heavy quark 
(here e c = 2/3). The quantities to, a s and / s stand for the bare mass, bare 
coupling constant and bare gluon density respectively. The meaning of the 
latter function will be discussed later. Further the order a™ contribution 

(n) 

to the heavy quark coefficient function will be denoted by H- k . Finally 
the superscript on the gluon density indicates that we are dealing with a 
singlet (S) quantity with respect to the flavour symmetry transformations. 
Here the flavour symmetry group is given by SU(n/). In the case of charm 
production we choose n/ = 3. 

The order a s corrections to the photon-gluon fusion process in Eq. (15) 
have been calculated in pf. Some of the virtual contributions to the Born 
reaction are shown in Fig. The corresponding Feynman integrals, denoted 
by n i ~ N {2tt)~ n J d N kf(k,ki), where k\ is an external momentum, reveal 
ultraviolet (UV), infrared (IR) and collinear (C) divergences. They arise 
when the integration momentum k takes the values k — ► 00 (UV), k = 
(IR) and k\\k\ (C) respectively. The former two divergences are very well 
known in quantum field theory and we refer the reader to the textbooks. 
The collinear singularities originate from the propagator l/(k — ki) 2 present 
in f(k, k\). It becomes singular when k 2 = k\ = so that {k — k\) 2 = —2 
I k || k\ I (1 — cos 9). The singularity then shows up at 9 = 0. The 
appearance of IR- and C-divergences is due to the fact that we neglect 
confinement effects in the initial and final state while doing perturbation 
theory. Here all external partons are put on-shell which is certainly not 
possible if we would have applied nonperturbative methods. In this case the 
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Fig. 5. The bremsstrahlungsprocess 7* +g — > c+c+g contributing to the coefficient 
functions H}„. 



partons are off-shell akin to the case when the particles constitute a bound 
state. Fortunately there exist some theorems in quantum field theory which 
enable us to remove the divergences mentioned above. This will be discussed 
below. For the moment we have to find a way to define the Feynman 
integrals in which these divergences occur. For that purpose we choose the 
method of iV-dimensional regularization which is the most suitable one since 
it preserves all Ward identities characteristic of gauge field theories. Using 
this method the divergences manifest themselves as pole terms of the type 
l/e k with e = JV-4. 

Besides the virtual gluon graphs in Fig. |4| there are also contributions from 
the gluon bremsstrahlung diagrams in Fig. §. This process is given by 

7* + g -> c + c + g. (17) 

Analogous to the one-loop Feynman integrals discussed above, IR and C di- 
vergences also show up in the phase space integrals corresponding to process 
(17). They arise when &2 — ► (IR) and ^H&i (C). The last condition orig- 
inates from the propagator f{k\, k%) ~ l/(&i — &2) 2 appearing in the phase 
space integral /u 4_7V J d N ~ l h.2{2E2)~ l /(&2 5 k\). Addition of the virtual cor- 
rections to process (15) and the contributions from the bremsstrahlung reac- 
tion in (17)) leads to a cancellation of all IR singularities. This is called the 
Bloch-Nordsieck theorem B. Furthermore the Kinoshita-Lee-Nauenberg 
theorem || states that collinear divergences which can be attributed to the 
final state cancel too. This cancellation only happens when the final state is 
completely inclusive. Notice that up to the order in a s discussed above final 
state collinear divergences are not present due to the heavy quark mass. The 
singularities which remain are of UV and initial state collinear origin. The 
former are removed by mass and coupling constant renormalization. Mass 
renormalization is performed by replacing the bare mass rh in the lowest 
order coefficient function (z,Q 2 ,m 2 ,a s ) via the substitution 
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Fig. 6. The Bethe-Heitler process j*+q^c + c + q contributing to the coefficient 
functions i?j„. 



'2 

rh 



m[l + a s (~6 + di)] with a s = ^-. (18) 

After the substitution one has to make an expansion around a s which in- 
volves taking derivatives of Hi )g with respect to m 2 . Further d\ is an arbi- 
trary constant which is fixed by choosing the on-mass-shell scheme. After 
this renormalization the second order coefficient function takes the following 
form. 



(2) _ - \r 1 1, ,m 2 ,,„ t ^ _ „m „ r 2 , ,m? 



{- + ~ 9 Inft)}PW » < - M— + ln(^)}< 
+tfg itc ' (2) . (19) 



Here we have distinguished between the UV and the C-divergences which 
are indicated by 1/euv and 1/ec respectively. Besides the higher order 
corrections to the photon-gluon fusion process there also exists another a 2 
subprocess where the photon couples to the charm quark. It is given by the 
Bethe-Heitler reaction (see Fig. ^). Instead of a gluon we have now one 
of the light (anti-)quarks in the initial state and the reaction proceeds as 
follows 

7* + q(q)^c + c + q(q). (20) 

In this process we only encounter collinear divergences and the bare heavy 
quark coefficient function reads 

h% = ^ + \ h^)}p$ ® 2$ + <r (2) ■ 

Corrected up to order a 2 the structure function reads 

Fi,c = e 2 c [P g ® {HQ + HQ } + fl ® HQ + -..]. (22) 
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Next we have to apply coupling constant renormalization in order to get rid 
of the remaining UV divergence. This is achieved by replacing 

a s = a s [l + a s p {^ + ln(^)+b 1 )] (23) 

in the coefficient function (z,Q 2 ,m 2 ,a s ). In the above expression [ir 
stands for the renormalization scale and a s = a s (ii 2 R ). Further (5q denotes 
the lowest order coefficient appearing in the beta-function which will be 
defined below. The constant b± is scheme dependent. Here we will choose 
the MS-scheme so that b± = je — ln(4-7r) (7^ is the Euler constant). Finally 
we have to perform mass factorization in order to get rid of the collinear 
divergence indicated by the pole term 1/ec in Eqs. (19), (21). This is 
achieved by replacing the bare fk by the renormalized parton density ft- In 
general one has 

ft = ftfrhA) ® [*« - a.P<?> {±- + \ ln(§) + cff)] , (24) 
where denotes the singlet combination of the light flavour densities 

rif 

/? = £[/* + /*]■ (25) 
k=i 

For rif = 3, the index k runs over (u, d, s). Further \ip stands for the 
factorization scale and pffi (k,l = q,g) denote the lowest order DGLAP 
splitting functions. The constants off are scheme dependent. Here we 
will choose the MS-scheme which implies that cm = Je — ln(47r). One can 
also perform mass factorization on the level of the heavy quark coefficient 
functions. In this case one gets the equation 



H i ,k(Q 2 ,m 2 ,a s ,fJ 2 R ,e c ,fJ 2 ) = 

r/fc(a s ,^>,/4,ec.,/i 2 ) ® H it i(Q 2 ,m 2 ,a s , n\, fi 2 F ) , (26) 
where denotes the transition function given by 

r« = 6 kl + a s P^ (1 + \ ln(§) + c«) + ■ ■ ■ . (27) 
From Eq. (24) we infer that 

fl(fi R ^ F )=r M (a s ,fi 2 R ,fi 2 F ,e c ,fi 2 )^fi i . (28) 
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Fig. 7. The Compton process 7* + q — > c + c + g contributing to the coefficient 
functions L,-„. 



In order to get finite coefficient functions up to order a 2 s we need the fol- 
lowing transition functions 

T qq = 1 + OK), 
T qg = 0(a s ) , 

1 1 - 2 
_ + o 



r M = l + a s Pi?(^ + ^ln(^) + c£)) + .-. (29) 

Notice that we have suppressed the z-dependence of the functions pf^' , 
and 1 = (5(1 — 2;). After renormalization and mass factorization the charm 
component of the structure function reads 

Fi, c = e 2 c [f 9 ® {H?J + H®} + f q ® H®] , (30) 

where H^g is obtained from H^g by the substitution : m — > m and a s — > a s . 
Further we introduce the shorthand notations 

< = <r (2) + a. (I H4)^ ® < - A< <£)) , (3D 



H% = <r >(2) + (5 ® O • (32) 

\z fj,p ' 

Besides the reactions discussed above where the photon interacts with the 
heavy quark the former can also couple to the light (anti-)quark. This 
happens for the first time in order o? s . In this order charm production 
proceeds via the Compton process (Fig. 0)) which is given by Eq. 20 
where now the photon couples to the light (anti-)quark. The corresponding 
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coefficient function is denoted by which in order a 2 is equal to L^g. 
The structure function Fj )C gets the contribution 

*i,c = ^"E e * [/I ® L L + n ffk S ® • (33) 



The non-singlet combination of the quark densities is defined by 

1 



fk S = fk + fk-—f? (k = u,d,s). (34) 



Up to second order no collinear divergences appear in this type of processes. 
However in higher order they do appear. Moreover there are reactions with 
a gluon in the initial state leading to the heavy quark coefficient function 
Lf g . Together with Lf they satisfy the same mass factorization relations 
as presented for Hf k in Eq. (26). For the non-singlet part a simpler relation 
holds (no mixing) which is given by 

fNS _ r NSrNS ( oe\ 
n i,q — 1 qq ^i,q ■ \°°) 

Likewise for the non-singlet parton density we have 

/f = r» s /f. (36) 

The reason that the transition functions are the same for and 
follows from the universality of collinear divergences. It means that the lat- 
ter are process independent. The same residues, represented by the DGLAP 
functions Pui , are also found while calculating corrections to other hard pro- 
cesses differing from heavy flavour production. One of the most important 
consequences is that the finite parton densities also become universal after 
mass factorization so that one can use them as input for other processes to 
yield absolute predictions. Collecting all contributions the charm compo- 
nent of the deep inelastic structure function for the proton reads 

1 n/ r 

F itC (n f ,Q 2 ,m 2 ) = — £ e\ /x 2 ) ® L% q (n h Q 2 , m 2 , /j 2 ) 

n f k=i 

+fg{n f , /x 2 ) ® Lfjrif, Q 2 ,m 2 , fi 2 ) 
+n / / fc NS (n / , /J 2 ) ® Lf*(n f , Q 2 , m 2 ,fi 2 ) 
+e 2 [fl (jif, [J> 2 ) ® Hfg (n f , Q 2 , m 2 , /j 2 ) 
+fS(n f , f j, 2 )®HtJn f ,Q 2 ,m 2 , f i 2 )} , (37) 
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Fig. 8. The function H 2 , g {x,Q 2 ,m 2 c ) for Q 2 = 10 (lower solid line), Q 2 = 100 
(lower dotted line), Q 2 = 10 3 (upper solid line), Q 2 = 10 4 (upper dotted line). All 
units are in (GeV/c) 2 



where now all quantities are finite. Notice that contrary to the functions 
^ are purely singlet (PS) only. This automatically holds for In the 
case of k = q we have added the superscript PS because = (no 

intrinsic charm !!!). As we have shown above the coefficient functions and 
the parton densities depend on two different scales (j,r and fip- In order 
to simplify the renormalization group equations (RGE's) below we will set 
them equal i.e. (Ir = (If = fJ-- Let us first define the total derivative with 
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Fig. 9. Same as fig. 8 for H 2 , q (x, Q 2 ,m 2 c ). 



respect to ji. 



d d n . . d 
D = fi— = n— + p{a s )—— , 
dji Ofi oa s 



(38) 



where (3(a s ) denotes the beta-function which can be written as a series 
expansion in the strong coupling constant 



(3{a s ) = -2p a 2 s - 2(3 ia 3 s + • • • . (39) 
One can show that the parton densities satisfy the following RGE's 

DC = P ? N >/f , Dfl = Pl^f^. (40) 
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Fig. 10. Same as fig. 8 for L 2 , q {x, Q 2 , ml). 

The heavy quark coefficient functions satisfy the RGE's 

UL i,q — ~^qq ® L i,q UL i,k ~ ~^lk L i,l 

DH* k = -i? fc ®fl?,. (41) 
The DGLAP splitting functions can be expanded in the coupling constant 
Pki = a s P^ +a 2 s P^ + ••• . (42) 

Using the above equations one can now show that Fj jC is a renormalization 
group invariant i.e. D Fi tC = 0. It means that it is a physical quantity 
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Fig. 11. Soft gluon contributions to the coefficient functions H ri 



(2) 



which is independent of the scale [i and the chosen scheme. (Notice that in 
these lectures we have adopted the MS-scheme for the renormalization of 
the coupling constant and for mass factorization). 

The order a 2 contributions to the heavy quark coefficient functions iTj j. and 
Li t k have been computed in Q. The expressions are so complicated that 
it is impossible to give analytic results. They are buried in long programs 
containing the numerical computation of two dimensional integrals. In order 
to make these coefficient functions more amenable for phenomenological 



applications they are tabulated [11] in the form of a two dimensional array 
in the variables 

ri=^~l , e=^- (43) 

The first variable is chosen in such a way that the threshold region r/ = 
(s — 4m 2 )(4m 2 )~ 1 ~ is exposed in a clearer way. Notice that this region 
dominates the integrand of the structure function (see fl2| ) 

F iyC (x,Q 2 ,m 2 )~ — f k (x/z)H Lk (z,Q 2 ,m 2 ). (44) 

Jx Z 



In Figs. g, □> |10| we have plotted the coefficient functions #2,31 -^2,9 and 
L>2,q respectively. Here we have chosen m = m c = 1.5 GeV/c. From these 
plots we infer that Hi,g is much larger than the two other heavy quark 
coefficient functions so that the photon-gluon fusion process constitutes the 
bulk of the radiative corrections. Since this process only depends on the 
gluon density fg(z,/j 2 ) charm electroproduction yields a measurement of 
the latter density, in particular at small z. The coefficient functions reveal 
large corrections which occur in two regions given by 

1 Threshold regime: s ~ 4m 2 or z ~ -z m ax = Q 2 /{Q 2 + 4m 2 ) 

Here we have two types of large corrections. The first one is due to 
soft gluon bremsstrahlung. As we have discussed before the real gluon 
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Fig. 12. The Coulomb singularity appearing in the box graph of the process 7* + 
c + c. 



and the virtual gluon corrections show IR divergences when &2 — > 0. In 
the real gluon case we have a phase space factor and the corresponding 
coefficient function behaves like 

ff S' REAL ~^(^) <m/1 <- («) 

On the other hand the virtual contribution equals 

<' VIKT ~ • («) 

em 

Addition of the two expressions above leaves the finite result 

H {2) ~lnr~ 4m2 V (1) (47) 

which however blows up when s — ► 4m 2 . The above behaviour is 
typical for QED. In QCD we have a C-divergence in addition to the 
IR divergence. Hence after having removed the IR and C-singularities 
the power of the logarithms will be doubled. The final result in QCD 
shows the following behaviour 

H {n) ~ a-ln 2 ^- 1 ) r~ 4m2 WW (48) 
The above logarithms can be summed in all orders of perturbation the- 



ory by exponentiating them (see [13|). The effect of these logarithms 



is shown in Fig. |3| where we plotted H^g- In the threshold region 
where x is large and Q 2 is small one observes a rise in the coefficient 
function. 

Another large correction near threshold can be traced back to the 



Coulomb singularity appearing in the virtual contribution in Fig. |12. 
One obtains the following expression 

| M (2),viRT |2 - -> #( 2 )< VIRT ~^ (49) 
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This singularity is enhanced by multiple gluon exchanges between the 
heavy quark lines. It arises because we apply perturbation theory. The 
correct way to deal with these exchanges is to apply nonperturbative 
methods. In this case the nonperturbative result shows a singularity at 
a s = which implies that one cannot make an expansion around this 
point. In practice the effect of the Coulomb singularity is so small that 
one can hardly see it in the heavy quark coefficient functions. item[2] 
Asymptotic regime: s 3> m? or z —* 0. 

The large corrections in this region are due to the exchange of multiple 
soft gluons in the t-channel. The corresponding coefficient functions 
have the form 

H^ n Mz,Q 2 ,m 2 ) ~ -ln n ~ 2 (z)h(Q 2 ,m 2 ) for n > 2; k = q,g. (50) 

For n = 2 see Fig. |l3[ Like in the previous cases these large corrections 
can be resummed (see [14|). The soft gluon exchange mechanism has 



a huge effect on the heavy quark coefficient functions. It is responsible 
for the large plateau in the small x region of the functions H 2 L (Fig. 

D and H^g (Fig. g). This plateau rises as Q increases. However it 
turns out that this production mechanism is not so important for the 
behaviour of Fi jC in Eq. (44) at small x since the the main contribution 
to the integral comes from large z rather than small z. 



4. Asymptotic heavy quark coefficient functions 

As we mentioned in the previous section it is very hard to get analytical 
expressions for the second order heavy quark coefficient functions except 



for Lj j(? (i = 2, L) which are published in [15|. However when Q 3> m 
it is possible to obtain analytical expressions. This is very useful because 
the latter serve as a check on the exact calculations carried out in p|, [|iT| . 
Furthermore it turns out that for charm production the asymptotic heavy 
quark coefficient functions give an equally good description as the exact ones 
when Q 2 > 20(GeV/c) 2 and x < 0.01. Finally the asymptotic expressions 
can be also used for the charm component of the structure function (Eq. 
37) in the so called variable flavour number scheme (VFNS) discussed in the 
next section. In order to compute the asymptotic heavy quark coefficient 
functions one can proceed in two ways. 

1 Standard method 



In this case one evaluates all Feynman and phase space integrals for 
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Fig. 13. Soft gluon exchanges in the t-channel of the processes: 7* + g 
and 7* + q — ► c + c + q. 



■c + j 



Q 2 S> ?n 2 . This method is still elaborate and one has to be careful 
with terms in these integrals which are proportional to (m 2 ) 3 because 
they survive in the limit indicated above. 

2 Techniques of asymptotic expansions 

Here we can use the operator product expansion (OPE) which is very 
elegant. As we will explain by presenting some examples one can 
compute the asymptotic expressions for Hi t k and Lj^ provided we 
know the massless parton coefficient functions denoted by Ci t k and 
the heavy quark operator matrix elements (OME's) given by A c k and 
Aki, c - (k,l = q,g)- Fortunately the functions Co- are already known 
up to order a 2 for some time [16|. Recently the OME's have also been 



computed up the same order (see [15], [17]) so that one can obtain the 
heavy quark coefficient functions without doing too much work. 

Adopting the second method in the subsequent part of this section it is very 
easy to show that one obtains the following result in lowest order 

<f MP '«(,, Q 2 ,m 2 ) = a s[P £\z) In (g) + h§>{z)] . (51) 

The above equation can be easily inferred from the literature [j|], ||] by 
taking the limit Q 2 3> m 2 of the exact formula. The latter can be also 
obtained in the limit m — > 0. In this case the charm mass acts as a regulator 
for the C-divergence which is indicated by its residue represented by the 
DGLAP splitting function P qg |l^]. Because of the analogy between lnm 2 
and the pole term 1/ec which represents the collinear divergence in the 
case of massless partons treated in the last section we can now apply mass 
factorization to remove them from the asymptotic heavy quark coefficient 
functions. In lowest order this is very simple and we obtain 



< SYMP ' (1) (.,Q 2 ,- 2 ) =C« ' (52) 
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with 



and 



C S0>^) =*.ir$H*)to$) +<£](*)], (53) 



4? (*, = a.[P& 0) (*) ^ (^) + ag>(*)] - (54) 

Here we want to stress the analogy between the quantities treated in the 
previous and the present section. They are 

section 3 section 4 



H i>k (z,Q 2 ,m 2 ,e c ) H^ YMP (z,Q 2 ,m 

Q 



Hi ik (z,Q 2 ,m 2 ,ii 2 ) C it k(z,-^) 

z, — _-). (55) 



In particular we want to emphasize the similar role played by the transition 
functions V k i and the OME's A ck . Both of them absorb the collinear diver- 
gences from the original quantities. We will now show that in lowest order 
the quantity A^j really stands for the heavy quark OME. 
Squaring the amplitude corresponding to the Feynman graphs in Fig. |l| 
and integrating over the final particle phase space one can apply the optical 
theorem which states that 

H it k(s,Q 2 ,m 2 ) = lmT itk (s,Q 2 ,m 2 ) , (56) 

where z = s/{s-\-Q 2 ). Here T^. stands for the forward Compton scattering 
amplitude (Fig |l4|) and the discontinuity denoted by Im is taken over the 
s-channel. If we define v = p • q = (s — m? + Q 2 )/2 one can write an 
unsubtracted dispersion relation which can be expressed as follows 

T iik (u,Q ,7/r) = - / _ dv — '— 



7T JQ2/2 V - V 

* dz' H itk (z',Q 2 ,m 2 



vr Jo z' 



1 00 rl 

-V2"" dz'z' n ~ l H it9 {z', Q 2 ,m 2 ) . (57) 
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Im k 



k + q 



k — p 

1 P IP* 
Fig. 14. Forward Compton scattering : 7* + g — > 7* + g and the optical theorem 



The graph for T^J in Fig. |l4| i yields the following integral 



r(i) 



^jT^(p)^(p) 

(2vr) w [(/c - p) 2 - m 2 }(k 2 - m 2 ) 2 [(k + q) 2 - m 2 



d kf(k,q,p,m). 



(58) 



where P/^ is a projection operator with i = 2,L and e^(p) stands for the 
polarization vector of the gluon. In the limit Q 2 3> m 2 the above expression 

ASYMP (1) 

will be called T- which can be split into the coefficient function 



-1 00 „i ^->2 

[ lim fi»kf(k, q ,p, m) - 

l Q 2 ^>m 2 J 



and the operator matrix element 



d k lim f{k,q,p,m) 

Q 2 ^>m 2 



(59) 



1 00 /"l 2 p 

-Y z~ n dz l z" n - l A^(z',^]= d N k lim /(*,<?, p, 

TT^ 7 9{ - VOL 2 ' J Q2 >>m 2 



m) . (60) 



The first term in [pASYMP,(i) ^ ag nQ singularity at m = and it carries the 



*,9 

whole Q 2 -dependence. The absence of the C-singularity can be checked by 
taking k oc p. The dependence on m is transferred to Aq^ which however is 
independent of q 2 = —Q 2 as we will show below. Further by interchanging 
limits and integrations in Eqs. (59) and (60) one introduces an artificial UV 

Notice that the lowest order Feynman graphs, shown in this section, are not complete. 
Current conservation and gauge invariance require additional graphs. 
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singularity which has to be subtracted before one gets the finite C^g and 
A^ig ■ To evaluate the expression in Eq. 60 we make the Taylor expansion 



(k + q) 2 -m 2 ~ q 2 ^ ,u \ -q 



= " E 7— 3vT >7mi7^m 2 • • • • • • , (61) 

n=l \ q > 



so that we get 



q ^ . . . q^n 



n=l n=l v ^ ' 

xif V(p)e» J ^ ^ )2 ) 1^2p?_ + m2 ]2 CT fe M2 ■ ■ • *W >(62) 

The vertex 2 n_1 7 Ml fe /12 • • • /c Mn originates from the heavy quark operator ma- 
trix element (OME) given by 

(c(k) \ O c ^..., n (0) \ c(k)) . (63) 

The heavy quark operator, which is a gauge invariant object (physical op- 
erator), is defined by 

O c , w ... Mn (x) = i>(x)-f^D^ ■ ■■D^{x) , (64) 

with the heavy quark (charm) field ij) and the covariant derivative = 

If we now sandwich the above operator between physical gluon states one 
obtains 

(g(k) | o w ...^(0) | 5 (fc)> = 



2 A / \ <T / \ / d N k 7A (^ +m)7 Ml (^ + m) 7(T 



,2 



~9fJ,ifJ,2 9fJ,3fJ-4P^5 ' ' ' Ppn ' ' ^symmetric • 

(65) 



Notice that the projection which survives in Eq. (62) is given by i = 2 
since contains the metric tensor g^ v only. Projections containing the 
vectors q^ and p^ vanish because of current conservation and the on-shell 
condition p 2 = respectively. Therefore j^ymp,(i) ^asymp,(i) 
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have no C-divergence for m = and the only contribution comes from C*jp 
which is Q 2 independent. The equality between Eq. (62) and Eq. (65) 

shows that the former represents the OME Acg as announced below Eq. 
(55). It is now easy to extract the OME if Eq. (65) is contracted by the 
tensor A w A M2 • • • A Mn where A^ is a lightlike vector i.e. A 2 = 0. Hence we 
get the expression 

A $ ,n d) = ii dz ' z,n ' ias [ p ^ {z,) ^ ■ (66) 

with 

2 i oo 2 

<G-4)=iE^'4r(^). (67) 

n=l 

After subtraction of the UV pole term, for which we choose the MS-scheme, 
one obtains the finite OME in Eq. (54). The finite coefficient function C^p 
Eq. (53) is obtained from Eq. (59) in a similar way. 
Generalizing mass factorization to higher orders one can write 



12, 

trASYMP^ 

tli 



L i,k 



\Q\ m *) = A lk (^)®C l , l (^ I ). (68) 



If we expand the OME's A^i and the light parton coefficient functions Cj & 
up to order a 2 , one obtains the following equations 

< SYMP ' (1) = 4><° c ) + 4°>< 1 ^ (69) 

C MP ' (2) = A $ ® Cg) + A® C&> + 45 cg> , (70) 

F ASYMP,(2) =4 2) 0C (O) +<)0C ( 2 ), 
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Here A gg is obtained by sandwiching the gauge invariant gluonic operator 
Og,fj,i—iMi\ x ) = F^D^z ■ ■ ■ D^ n _ 1 F tln ^ a , (72) 

between physical gluon states The heavy quark coefficient functions La- 
have to be added to C^fc before we can apply mass factorization 

M n /'fJ) +^f YMP (Q 2 V) = 4*, c (^) ®Cy(n/ + l,^) • (73) 

This relation also involves a redefinition of the strong coupling constant 
which in the right-hand-side now depends onnj + 1 instead of rif. The main 
effect of this operation is that the number of light flavours nj appearing in 
Cik is enhanced by one unit. The Aki yC are given by all light parton OME's 
containing charm loop contributions to the gluon self energy only. Up to 
order a 2 we have 

L ASYMP,( 2 ) = ^ C (0) + A (0) 3 C ( 2 ) + as/?0jc C W ln X ) . (74) 

A* 

Finally we want to emphasize that beyond the second order the above mass 
factorization relations become much more complicated. Here the heavy 
quark OME's A c ^ can also contribute to Eq. (73) whereas the OME's c 
also enter Eq. (68). Using the mass factorization relations in Eqs. (69)- 
(71) we have computed the asymptotic heavy quark coefficient functions 
H^ YMF '^ for k = q,g and i = 2,L in |l5|]. The same was done for 

L ASYMP,(2) (?4) _ Thig 

was possible because the light parton coefficient 



functions and the heavy quark OME's were already known in [16] and [15] 



respectively. The expressions reveal the following characteristic behaviour 



Ht^ {z ,Q\ m \f)~a{ £ QtfQ, (75) 

1 1 L ' K lib ' 

n+j<( 



with an analogous expression for L^? YMP . One of the most important fea- 



tures, discovered in [15|, is that the asymptotic expressions tend to the exact 
heavy quark coefficient functions at rather low £- values (see Eq. (43)) pro- 
vided z is not too small. According to ref. |0| £ = 10 for z < 0.01 which in 
charm production corresponds to Q 2 = 22.5 (GeV/c) 2 with m = 1.5 GeV/c 
The consequences of this behaviour for the charm component of the struc- 
ture function will be discussed in the next section. 
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Fig. 16. Two-loop vertex correction containing a heavy quark loop. It contributes 
to d 



5. Description of Fj iC in the various schemes 

In the literature one has proposed various descriptions of the charm 
component of the structure function in the framework of extrinsic charm 
production. Here one can distinguish the following schemes. 

1 The three flavour number scheme (TFNS) 

Here the production mechanisms are given by the photon-gluon fusion 
process and the higher order reactions. The charm component of the 
structure function Fi >c is given by Eq. (37). The parton densities are 
given by the three light flavour densities u,d,s and the gluon density 
9- 

2 The four flavour number scheme (FFNS) 

In this case Fj jC is expressed into convolutions of light parton coefficient 
functions with light parton densities. The latter are represented by four 
light flavours, which includes the charm quark, and the gluon. 

3 The variable-flavour number scheme (VFNS) 

This scheme interpolates between the results of the structure functions 
F ijC obtained from the TFNS and the FFNS. 

In this section we will now derive the expressions for i^ )C in the FFNS and 
the VFNS from the TFNS . 

In order to derive Ff™ s we need the deep inelastic structure function which 
receives contributions from subprocesses containing light partons in the ini- 
tial and final state only. It is given by 



Fi(n f ,Q 2 ) = 
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^ E 4 [/?(»/ , /x 2 ) ® (n,, ^) + /f (n /} Cf^ (n f , ^) 



-n f f™(n f ,^)®C?S(n f ,^)]. (76) 



In the case of TFNS we have to choose nj = 3. In the above expression 
the coefficient functions contain heavy quark loop contributions to the 
gluon self energies appearing in the light partonic matrix elements. This 
happens for the first time in order a 2 . An example is shown in Fig. 16. 
The charm quark loop contributions will now be removed from Fi(nf,Q^) 
and added to Fi c in (37). In the subsequent part of this section the latter 
will be called F^ C XACT . In this way the coefficient functions Li k get the 
proper asymptotic behaviour when Q 2 3> m 2 as indicated in Eq. (75). This 
allows us to perform mass factorization after we have added the light parton 
coefficient function Cj^ in which all contributions from heavy quark loops 
are removed (see Eq. (73)). Next we define the following quantity 



F l f YM "(n f ,x,Q 2 ,m 2 )= lim Fj c (nf,x,Q,m) . (77) 

Notice that i^ ASYMP j g g^ ven by the same expression as -Fj E c XACT except 
that now the exact heavy quark coefficient functions are replaced by their 
asymptotic analogues which have the form presented in Eq. (75). In Fig. 



17 we have plotted in next-to-leading order (NLO) the ratio 



F 2 ASYMP (x,Q 



2 m 2 



fi2 (.,QW) = - A ^-— . (78) 

For this plot we have adopted the parton density set with A4 = 200 MeV 
from ]ll| (MS-scheme). The mass factorization scale is chosen to be /i = Q. 
The charm quark mass is equal to m c = 1.5 GeV/c. From this figure we infer 
that for Q 2 > 20 (GeV/c) 2 and x < 0.01, F 2 A C SYMP coincides with F 2 EXACT . 
For Fl, c this happens when Q 2 is much larger i.e. Q 2 > 1000 (GeV/c) 2 . 
This implies that the large logarithmic terms in the heavy quark coeffi- 
cient functions H^k and L^k as given by Eq. (75) entirely determine the 
charm component of the structure function. Since these corrections vitiate 
the perturbation series when Q 2 gets large they should be resummed in 
all orders of perturbation theory. This procedure has been carried out in 
|fl7| and it consists of four steps. First we add Fj(n/,Q 2 ) in Eq. (76) to 
Ff^ YMP (n f , x , Q 2 , m 2 ) in Eq. (77) and choose rif = 3. Second we apply 
mass factorization to the asymptotic heavy quark coefficient functions ac- 
cording to Eqs. (68), (73). In the third step we define new parton densities 
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Fig. 17. i?2(NLO) plotted as a function of Q 2 at fixed x; x = 0.1 (dashed-dotted 
line), x = 0.01 (dotted line), x — 10~ 3 (dashed line) , x = 10~ 4 (solid line). 



which are now presented in a four flavour number scheme (FFNS). For the 
three light flavour densities k = u, d, s we have 



2 

/ fc (4,/. 2 ) + /s(4,/i 2 ) = 4^ s c (3,^) ® [/ fe (3,M 2 ) + /fe(3, 

+^c(3,^)®S(3 > m 2 ) 
2 

+^( 3 '^)® G ( 3 .M 2 )- 
The gluon density in the FFNS reads 

G(4, /. 2 ) = ^ g c (3, ,u 2 ) ® E(3, ^ 2 ) + 4 PiC (3, m 2 ) ® G(3, » 



(79) 



(80) 



Finally we get a new light flavour density represented by the charm quark 
/ c+c -(4,/i 2 ) ee / 4 (4, / u 2 ) + / I (4, / u 2 ) 

= ^(3,^)0S(3, M 2 )+^(3,^)0G(3,/, 2 ). (81) 
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Fig. 18. Scale dependence of the structure functions at x = 0.05 and Q = 10 GeV/c; 



F, 



VFNS,(2) 



(solid line), F 2 P ° P ' (2) (dashed line), ^ ACT '( 2 ) (dotted line). 



The above charm quark density has the property that it does not vanish 
at fi = m in the MS-scheme contrary to what is usually assumed in the 
literature (see e.g. @, @, 0]), ||). In the fourth step we rearrange 
terms and obtain 



F- 



ASYMP 



i,c 



(n f ,x,Q 2 ,m 2 )+F l (n f ,x,Q 2 ) = F i (n f + l,x,Q 2 ), (82) 



which is the FFNS result for the total structure function. From the latter 
quantity one can extract the expression for the charm quark component of 
the proton structure function in the FFNS which will be denoted by 



i£ DP (n / + l,Q 2 ) 



f c+5 (n f + l,iS)®Cf!°(n f + 1, 2_J + E(n, + 1,//) 8> 



CfS(n f + 1, ^) + G( n/ + 1, ^ 2 ) ® Cf, 9 ( 



/i- 

where we have defined 

Cf, q (n f )=Cf^(n f )+n f C^(n f ) 



1,9 



n f C l 



(83) 



(84) 
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The superscript PDF in Eq. (3) stands for parton density function which 
means that the charm component of the structure function is completely ex- 
pressed into parton densities multiplied by the light parton coefficient func- 
tions. Notice that Ff^ F is a renormalization group invariant like Ff^ T 
and F/^, SYMP so that they satisfy the equation D Fj jC = (see Eq. (38) and 
below). Further -F P C DF originates from the charm quark coefficient functions 
Hi & from which it follows that the former is proportional to e 2 only. The 
functions Ljfc, which are multiplied by the light charge squared e|, con- 
tribute to both FfP F (nf + 1, Q 2 ) and F 2 (nj + 1, Q 2 ) so that the number of 
flavours in these structure functions are increased by one unit. The FFNS 
charm quark density is mainly determined by the size of the TFNS gluon 
density G(rif,z,n 2 ) for nj = 3. Therefore the latter plays a major role in 
the behaviour of Ff^ as well as of i 7 f c DF • An analysis of both structure 



functions in [17] reveals that the former gives the best description of charm 
production in the threshold region where Q 2 is small and x is large. On 
the other hand when Q 2 is large and x is small it turns out that it is better 
to use F^ F because it is in this region where the large logarithms in Eq. 
(75) dominate so that they have to be resummed. Therefore the TFNS is 
the most suitable scheme for the charm component of the structure function 
near threshold whereas far away from this region it turns out that the FFNS 
is more appropriate. 

One also needs a scheme which merges the advantages of these two 
pictures and provides us with good description of F 2iC in the intermediate 
regime in Q 2 . This is given by the so called variable flavour number scheme 
(VFNS). In |l7| we proposed the following VFNS structure function 



F^(x,Q 2 ,m 2 ) = F^ F (n f + l,x,Q 2 ) + F^(n f ,x,Q 2 ,m 2 ) 

-F^™ p (n f ,x,Q 2 ,m 2 ). (85) 

The above expression is a generalization of Eq. (9) in ||23[|, which was only 
presented in leading order (LO) and has been implemented in a recent global 
parton density analysis p^| . (A different VFNS scheme has recently been 
proposed in [24|.) In LO the VFNS scheme has the properties that for 



Q 2 » m 2 , F™ CT -> F 2 A C SYMP which means that F 2 VFNS -» F P C DF Further 
at low Q 2 (i.e. Q 2 < m 2 ) F 2 ASYMP - F P C DF so that F 2 VFNS F 2 FXACT . 
provided we put z max = 1 in i ? 2 A ( ? YMP . However the last relation is no 
longer true in higher order in ot s . The main reason is that new production 
mechanisms appear in NLO giving rise to the coefficient functions L^. The 
latter contain the prefactor e| (e.g. the Compton process) and show up in 
F 2 EXACT and F^ YMP but not in F F c F , which is proportional to e 2 c only 
(see Eq. (83). This higher order effect is only noticeable in the threshold 
regime where x is very large and Q 2 is very small. 
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Fig. 19. Structure functions at Q 2 = 25 (GcV/c) 2 ; c "°'^> (solid line), F^ c 
(dashed line), ^XACT.p) ^ ^ e( j i me ) T^he experimental data are from [6] (closed 
circles) and [7] (open circles). 



As an application we have studied the charm component of the proton 
structure function in the three schemes mentioned above. The coefficient 
functions used for these structure functions are all computed up to order a 2 

(see f§, (T^]) so that we will denote them by F 2 E f ACT ' (2) (Eq. (37)), 
f pdf,(2) ^ ^ 83 ^ and f vfns,(2) ^ i nFig |JJ we have plotted the 

scale (//) dependence up to NLO of the structure functions in the different 
schemes mentioned above. For these plots we used the parton densities 
presented in the MS-scheme with A4 = 200 MeV in [20|. The scale fi is 
adopted from pjj and it is given by 



yf= m 2 + kQ 2 (l 
- 2 for 



m 



m 2 /Q 
Q 2 < m 2 



2\n 



for 



Q 2 > m 2 



(86) 



with k = 0.5, n = 2 and m = 1.5 (GeV/c 2 ). In [§3| and in g| it was 
shown that in LO is less sensitive to variations in the scale \x than 

each term on the right-hand-side of Eq. (85) separately. However in NLO 
we observe in Fig. |l^ that there is no reason to prefer one scheme over 
the other. This is corroborated by the findings in [^FJ] and 12] where one 
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x = 0.0005(X4000) 




10 Q 2 (GeV 2 ) 



Fig. 20. Structure functions ^fns,(2) jj nc ^ i^° F '( 2 ) (dashed line), 

j^xact,^) (dotted line). The experimental data are from [27], [6] and [7]. 



EXACT (2) 

could show that there is a considerable improvement in F 2c ' with 
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respect to variations in the mass factorization scale when this quantity is 



computed up to NLO. In Fig. |l9|we plot i*2, c hi the three different schemes 
and compare the results with the recent data from the HI || and ZEUS @ 
collaborations. Further we have chosen the parton density set in |l9| with 
A 4 = 200 MeV. Notice that all parton density sets in j20|, used 
for the plots of the structure functions are presented in FFNS although in 
principle one has to choose a TFNS parametrization for both the parton 
densities (see e.g |l8|) and the running coupling constant for the compu- 



tation of F^^^ and F^ YMP . However we have checked that the latter 



are not significantly altered when we replace the parton densities in |18| by 
those in JlEf] . From Fig. 19 we infer that the data are in agreement with 



all schemes in which the structure functions are computed. Further the re- 
sults for i? v J NS '( 2 ) are always between j?exact,(2) anc [ F^ F '^ 2 \ it is clear 

that one needs more precise data in finer bins of x and Q 2 to discriminate 
between the various schemes. Finally we have made a comparison with the 
data in the small x-region as well as in the large x-region obtained by the 
experimental groups HI |], ZEUS @ and EMC §7) respectively. The plots 



for the various schemes are shown in Fig. 2C. The agreement between the 
data and the predictions from the structure functions in Eq. (85) is fairly 
good except for x = 0.237. The discrepancy occuring in the threshold region 



is mainly due to the large negative contribution coming from Fig. 16 and 



the chosen scale in Eq. (86) which originates from 22]. It turns out that at 
low Q 2 this scale becomes too small so that the running coupling constant 
is too big. In this region all schemes lead to negative structure functions in 
particular for CT . Therefore perturbation theory breaks down and one 
has to choose a larger scale. This phenomenon, which was not observed in 
the LO analysis in pSfl , only appears when the NLO corrections are taken 
into account. 
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